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Abstract 

We consider the Poisson-Boltzmann equation in a periodic cell, representative of a porous medium. It 
is a model for the electrostatic distribution of N chemical species diluted in a liquid at rest, occupying the 
pore space with charged solid boundaries. We study the asymptotic behavior of its solution depending on a 
parameter /3 which is the square of the ratio between a characteristic pore length and the Debye length. For 
small /3 we identify the limit problem which is still a nonlinear Poisson equation involving only one species 
with maximal valence, opposite to the average of the given surface charge density. This result justifies the 
Donnan effect, observing that the ions for which the charge is the one of the solid phase are expelled from the 
pores. For large (3 we prove that the solution behaves like a boundary layer near the pore walls and is constant 
far away in the bulk. Our analysis is valid for Neumann boundary conditions (namely for imposed surface 
charge densities) and establishes rigorously that solid interfaces are uncoupled from the bulk fiuid, so that the 
simplified additive theories, such as the one of the popular Derjaguin, Landau, Verwey and Overbeek (DLVO) 
approach, can be used. We show that the asymptotic behavior is completely different in the case of Dirichlet 
boundary conditions (namely for imposed surface potential). 
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1 Introduction 



Originally proposed at the dawn of the XX*'^ century by Gouy and Chapman [U [5], the Poisson-Boltzmann 
(PB) equation is still the corner stone of most of the theoretical descriptions of electrokinetic phenomena. Many 
works emphasized the limitation of such a model in the last decades, though. The ions are only represented by 
their charge, they do not have any volume, the correlation are neglected. The molecular nature of the solvent 
and further specific forces (such as the London dispersion) are completely ignored [H [9]. Thus the domain of 
validity appears to be relatively narrow, typically in the regime of dilute simple (most of the time monovalent) 
electrolytes. Nevertheless, because of its simplicity, most of the theories of equilibrium and transport in charged 
diphasic media are still direct generalization of the PB approach. For example, geological media (such as clays) 
[151 [H], electrochemistry [3l[10], and colloidal physics [13] are still based on the original concepts described by 
the Poisson-Boltzmann equation. 

The success of such an approach is due to several aspects. It justifies the popular Derjaguin, Landau, Verwey 
and Overbeek (DLVO) theory j21j that explains the stability of charged suspensions. In the case of charged 
porous media, the Poisson-Boltzmann approach is also particularly significant because it yields the equilibrium 
electrostatic properties of the materials and it can be easily coupled to further equations in order to provide 
a global model of the system. Indeed, for the transport properties, the Poisson-Boltzmann equation can be 
extended in order to give the Poisson-Nernst-Planck (PNP) formalism which describes non-equilibrium processes 
in complex systems [23]. For example, in the case of clays, the description of electrokinetic processes in the large 
pores (meso and macroporosities) can be performed thanks to the PB equation [T7]. The molecular nature of 
the system are found to be important only for micropores (typically for distances less than 2 nm) |14] . 

In porous media, the PB exhibit two different regimes, depending on the value of the salt concentration or 
the pore sizes on the system. 

• If the pore size L is much larger than the Debye length of the electrolyte, the solid charge is screened 
by the microscopic ions. Thus, the local charge density is globally zero, but at the interface. Because 
of the relatively small value \o, this case corresponds to numerous applications. The solid interfaces are 
uncoupled so that the DLVO approach is valid. Far away from the interface, the coulombic forces can be 
modelled by effective parameters, such as the effective charge [20] or the zeta potential. 

• Conversely, if the pore size L is much smaller than Xo, the charge of the solid surface is not screened. 
It means that the resulting electrostatic force is important anywhere in the material. Ions for which the 
charge is the same than the one of the solid phase are expelled from the material (Donnan effect) [9]. 
The electro-osmotic flux becomes especially important. This case is significant because it corresponds to 
nanoporosities at low salt concentration. 

The two asymptotic limits can be taken into account thanks to the coupling parameter (5 = (L/Xd)"^. The large 
pore size domain (large /3) corresponds to most of the porous systems. Nevertheless, many porous materials 
exhibit microscopic pores for which the opposite limit (small /3) is relevant. For example, montmorillonite clays 
have different porosities, and the smallest ones, which are obtained at very low hydration, are even of the order 
of molecular distances. 

The goal of the present paper is to give a rigorous mathematical analysis of these two opposite asymptotic 
limits. The paper is organized as follows. Section [2] introduces the model and defines the relevant reduced 
units. Section [3] describes precisely the geometry of the porous cell and discusses the issue of the existence and 
uniqueness of the solution to PB equation. Section U] studies the limit case of very small pores, i.e., when f3 goes 
to zero. Section [5] is concerned with the opposite situation of very large pores, when /3 goes to infinity. Eventually 
Section [6] investigates the case of Dirichlet boundary conditions (namely for imposed surface potential) instead of 
Neumann boundary condition (namely imposed surface charge). A brief description of our main results is given 
in the next section after introducing the necessary notations. 

2 The model and our main results 

We consider the Poisson-Boltzmann system which describes the electrostatic distribution of N chemical species 
diluted in a liquid at rest, occupying a porous medium with charged solid boundaries. The electrostatic potential 
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is calculated from the Poisson equation 

N 

£A^* = -e^ZjU* in the bulk, (1) 
i=i 

where £ = £o£r is the dielectric constant of the solvent, e is the electron charge and 1 < j < iV, are the 
species concentrations. Since the pore walls are charged, the corresponding boundary condition is of Neumann 
type 

£V"^* ■ n = —a* on the surface, (2) 

where a* is a given surface charge and n is the unit exterior normal. 

At equilibrium the species concentrations n* are given by the Boltzmann distribution which corresponds to 
a balance between the chemical potential and the electrical field 

V{kBTlnn*) = -V{zje^*). (3) 

where Zj is the valence of the j-th species, the Boltzmann constant and T the temperature. It follows from 
([3]) that there exist positive constants n*(oo) (called infinite dilution concentrations) such that 

n* = n*(oo)exp|-^|. (4) 

The Poisson-Boltzmann system is the combination of ([1]) and (jH), together with the boundary condition ([2]). 

In order to make an asymptotic analysis of the Poisson-Boltzmann system, we first adimensionalize equations 
©5 ©1 &■ We denote by L the characteristic pore size and by ric the characteristic concentration. We introduce 
the Debye length defined by 



Ad 

and define a parameter 



£kBT 

2 ' 



/? = i^f. (5) 

The parameter /3 is the fundamental physical characteristic which drives the transport properties of an electrolyte 
solution in a porous media. For large /5 the electrical potential is concentrated in a diffuse layer next to the 
liquid/solid interface. Co- ions, for which the charge is the same as the one of the solid phase are able to go 
everywhere in the pores because the repelling electrostatic force of the solid phase is screened by the counterions. 
The electrostatic phenomena are mainly surfacic, and the interfaces are globally independent. For small /3, co- 
ions do not have access to the very small pores (Donnan effect). The local electroneutrality condition is not valid 
anymore and the electric fields of the solid interfaces are coupled. 

Furthermore, we define a characteristic surface charge density ac by 

£kBT 
eL 

and adimensionalized quantities 

a* ^ e"^* n* n*(oo) 

CJc kBT ' Uc ■' Uc 

Rescaling the space variable y = x/L, this yields the adimensionalized Poisson-Boltzmann system 

N 

Ay* = -/3^Zjnj(*) with nj(^') = n°(oo)e"^J* in the bulk, (6) 

and 

Vy* • n = —a on the surface. (7) 
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The goal of the present paper is to study the hmit of equations ^ and ([7]) when the parameter f3 is either very 
small or very large. Section 13.21 gives a precise mathematical framework for the Poisson-Boltzmann system. To 
discuss our results we sort the valences by increasing order and we assume that there are both anions and cations 

Zl < Z2 < ■■■ < zj\r and zi < < zn- 

Section m is devoted to the asymptotic analysis of ([6])- ([7]) when f3 goes to zero. This case corresponds to very 
small pores, L « Xd. In view of the definition of the Debye length A/j, a small value of f3 corresponds also 
to a small characteristic concentration ric- The asymptotic regime depends on the sign of the averaged charge 
Jg a dS. If it is negative (which means that the surface is positively charged), then Theorem 1111 states that only 
the anion with the most negative valence [zi] is important and that the potential behaves as 

^ ~ h v^o, 

where (pQ is the solution of the reduced system, involving only the species 1, 

r A(/?o = -^in?(cx))e-^i'^« in the bulk, 
\ V(/5o • n = —a on the surface. 

Note that the constant \og(i/z\ is going to +oo since zi < 0: it is a manifestation of the singularly perturbed 
character of the asymptotic analysis. As a consequence, the cation concentrations goes to zero while the ion 
concentrations blow up as rij = 0(/3~^-'/^i) and ni >> rij for j ^ 1- Of course, a symmetric behavior (involving 
only the species with most positive valence z^) holds true when the averaged charge is positive. On the other 
hand, when the averaged charge is zero, then the limit problem is much simpler (there is no singular perturbation) 
and it is given by Theorem [TH 

Section [5] is concerned with the opposite situation when /3 goes to infinity. This case corresponds to very large 
pores, L >> A/) or to large values of the characteristic concentration Uc- Such a situation is well understood in 
the physical and mathematical literature: it gives rise to a boundary layer (the so-called Debye's layer) close to 
the surface. For example, it has been analyzed in [2], [16] but the analysis is restricted to 1-d or similar simplified 
geometry. Our main contribution. Theorem [20| derives this boundary layer in a general geometric setting and 
gives a rigorous error estimate. Locally, close to the surface, the potential behaves as 



^(y) w , = exp 

^JPEUz]nOioo) 

where d{y) is the distance between the point y and the surface. Away from the surface, the concentrations nj 
are constant and satisfy the so-called bulk electroneutrality condition. Some numerical simulations for varying 
/3 can be found in our other work [1]. 

In the literature, one may find instead of the Neumann boundary condition ([2]) a Dirichlet one 

\I' = C on the surface, (8) 

where C is an imposed potential. Eventually Section [6] investigates the case when the Neumann boundary 
condition ([2]) is replaced by the Dirichlet boundary condition ([8]). Lemma [22] gives the asymptotic behavior for 
small /3: the analysis is quite obvious since there is no more singular perturbation. The potential now behaves 
as 

where ^'i is the solution of the corrector problem 

( ^ 

A^i = - ^ Zjn°(oo)e~^j '^ in the bulk, 

i=i 

^1 = on the surface. 









/3^z|nO(oo)| 
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Theorem 1281 gives the asymptotic behavior for large (3. It is again a boundary layer but with a totally different 
profile. More precisely we establish 

M/(y)«^o,c(v^%)) 

where d{y) is the distance between the point y and the surface and ^'o,( is the solution of the nonlinear ordinary 
differential equation (|113p . solution which, starting from the boundary value (" on the surface, is exponentially 
decaying at infinity. 



3 Geometry and existence theory for the Poisson-Boltzmann equation 

In this section we define the geometry and recall a classical result on the existence and uniqueness of the solution 
of the Poisson-Boltzmann equation. Furthermore we establish L°°-bounds for the solution. 



3.1 Geometry of the porous cell 

To simplify the presentation we consider a sample of a porous medium which occupies the periodic unit cell 
Y = (0, 1)'^ (identified with the unit torus T'^). The space dimension is typically d = 2,3. The pore space Yp, 
containing the electrolyte, is a 1-periodic smooth connected open subset of Y. More precisely, we consider a 
smooth partition Y = Ys UYp where I5 is the solid part and Yp is the fiuid part. The liquid/solid interface is 
S = dYs\dY. It is known that for a general Yp with non-empty boundary S, the distance function d, defined by 
d{y) = dist {y, S), y G Yp, is uniformly Lipschitz continuous. Let Yp = {y & Yp : d{y) < /Lt}. If we assume S to 
be smooth of class C^, then, for sufficiently small /x, the distance function has the same regularity, d G C^(Y^). 
Furthermore, by the Collar Neighborhood Theorem (see e.g. [H]), there exists a tubular neighborhoods U, 
S C U, isomorphic to 5 x /, for / C M an open interval. 

For yo G S, let n(?/o) and T(yo) denote respectively the unit normal to S at yo (exterior to Yp) and the tangent 
hyperplane to S at yo- Without loosing generality, we can suppose that in some neighborhood A/" = Af{yo) of yo, 
S is then given by yd = 7(y')> where y' = {yi, . . . ,yd-i), 7 G C^(T(yo) n AA) and ^y'^iy'o) = 0. The unit normal 
on S, at the point y = {y',^{y')) is 

Vl + |V».7(!/')P 

Here Ys corresponds to yd > 7(y')- 

Following jTj, page 355, for each point y £ Yp, there exists a unique point z = z{y) G S such that \y—z\ = d{y). 
The points y and its closest point z £ S are related by 

y = z — n{z)d. (9) 

By the Inverse Mapping Theorem applied in a neighborhood of yo £ S {S can be covered with a finite number 
of such neighbourhoods), dH) defines uniquely a principal coordinate system {q',qd) with q' = z' and qd = d{y) 
which are functions of y. The first coordinates q' are the tangential coordinates to S while the last coordinate 
qd is the signed distance to S taken with the sign -1- in the interior of Yp and with — in its exterior. We denote 
by = the local representation of the boundary S in the neighborhood of y^. Of course, the distance function 
d £ satisfies Vd{y) = —n{q'{y)) and \\/d{y)\ = 1. Furthermore, the tangential coordinates q' = {qi, . . . ,qd-i) 
can be chosen in such a way that the Hessian matrix [Z)^7((/'(?/o))] is diagonal at yo. For more details we refer 
to [7], pages 354-356. 



3.2 Mathematical framework for the Poisson-Boltzmann equation 

We start from the adimensionalized Poisson-Boltzmann system ([6]) and ^ (where we drop the index y for the 
spatial derivatives) that we complement with periodic boundary conditions on the outer boundary of the unit 
cell Y. More precisely, the potential ^ is a solution of the Poisson-Boltzmann equation 

-A'^ + fS^i^) = in Yp, 

V* • n = -0- on S, (10) 
^ is 1 — periodic, 
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where cr{y) is the adimensionaUzed charge distribution on the pore surface, /3 > is a parameter defined by 
([5]) as the square of the ratio of the pore size and the Debye length, and <I>(^') is the bulk density of charges, a 
nonlinear function defined by 

TV 

$W = -J^z,n,W (11) 
i=i 

where rij is the concentration of species j given, at equilibrium, as a function of the potential 

nj{^) = nO(oo)e-"^* (12) 

with a constant n^(oo) > 0. We assume that all valencies Zj are different. If not, we lump together different ions 
with the same valence. Of course, for physical reasons, all valencies zj are integers. We rank them by increasing 
order, 

Zl < Z2 < ... < ZN, 

and we denote by and j~ the sets of positive and negative valences. We also assume that they are both anions 
and cations, namely positive and negative valences, i.e., < < zat. The function $ is monotone since it is the 
derivative of a convex function 

N 

$(*) = C'(^) with C(^) = ^n°(oo)e-"^*. (13) 

The boundary value problem (jlOp is equivalent to the following minimization problem: 

inf Sj{^) = l[ \Vv\^dy + (3[ C{(t>)dy+ [ a^ds], (14) 

•peH^^iYp) I 2 Jyp Js J 

with H}^{Yf) = {93 G H^iYp), is 1 — periodic}. The functional J is strictly convex, which gives the uniqueness 
of the minimizer. Nevertheless, for arbitrary non-negative /3, n^, J may be not coercive on H^iYp) if all Zj's have 
the same sign. Therefore, we must put a condition on the Zj's so that the minimization problem (jl4p admits a 
solution. Following the literature, we impose the bulk electroneutrality condition 

N 

_$(0) = J^z,n°(oo) = 0, (15) 
i=i 

which guarantees that for o" = 0, the unique solution of (|10p is ^ = 0. Under assumption (|15p it is easy to see 
that J is coercive on H}^{Yp). 

Remark 1. The hulk electroneutrality condition M5\l is not a restriction. Actually all our results hold under the 
much weaker assumption that all valences Zj do not have the same sign. Indeed, if i fi5|) is not satisfied, we can 
make a change of variables in the Poisson-Boltzmann equation defining a new potential ^' = ^' + where 

is a constant reference potential. Since the function <I>(^''^) is continuous and admits the following limits at 
infinity 

lim ^>(^°) = ±00, 

there exists at least one value such that (^{^^) = 0. This change of variables for the potential leaves IT^J 
invariant if we change the constants n^(oo) in new constants h^{oo) = nj(oo)e~^J * . These new constants satisfy 
the bulk electroneutrality condition il5\) . 

In the sequel we assume that the charge density is a continuous periodic function 

a{y) G C#{S). (16) 
The well-posedness of the boundary value problem is a classical result. 
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Lemma 2 (|T2]). Assume il6\} and that the electroneutrality condition i fi5|) holds true. Then problem {1^^ has 
a unique solution ^ G H^lYp) such that 



N N 



^Zj-e"^^* and^Y^ 



zje 



are absolutely integrable. 



The issue of the solution boundedness was not correctly addressed in [12] . where it was merely proved that 
^ € U'{Yp) for every finite p. By using elementary comparison arguments we will prove a L°°-estimate in 
Proposition [3] (a similar result is also proved in [6]). To this end we introduce the following auxiliary Neumann 
problem 

-AC/ = [ a dS in Yp, 

XJTT ' ' (17) 

V(7 • n = —a on 6, ^ ' 

[/ is 1 — periodic, Jy^ U{y) dy = 0. 

Remark that (I17p admits a solution U € H^iYp) since the bulk and surface source terms are in equilibrium. 
Furthermore, the zero average condition of the solution gives its uniqueness. Under condition ()16p it is known 
that U is continuous and achieves its minimum and maximum in Yp- 
Then our L°°-bound reads as follows 



Proposition 3. The solution ^' of ( fjOj) satisfies the following bounds 

1 la -s-^ ZjU^Aoo) \ 

Uiy) - U„, - - logma. (^l. ^^-^ - ^ J > > 

1 la fr^ ZjU^Aoo) \ 

Uiy) - Um - — logmax 1, - J] ' ^ , (18) 

ZN I PzNn'}^{(X)) ZNn)^{oo) j 

where the sets and j~ denote the sets of positive and negative valences, respectively, and 

a = — — j- I a dS , Um = min Ujy) and Um = max U (y) . 
\Yf\ Js yeYp yeYp 

Proof. We use the variational formulation for ^ — U, which reads, for any smooth 1-periodic function ip, 

/ V(* - [/) • V(/3 + /? / ^{^)ipdy + a if dy = 0. (19) 

Jyp JYp JYp 

We take '^{y) = {"^{y) — U{y) + C)~ , where C is a constant to be determined and, as usual, the function 
/~ = min(/, 0) is the negative part of /. We note that by Lemma [21 ip G H^{Yp) and ^{'^)ip is integrable. It 
yields 

/ iVifl"^ dy + 13 [ {'^{^f)-<^{U-C)){^-U + Cydy 

JYp JYp 

+ [ i/3<l>iU -C) + a)pdy = 0. 

JYp 

By monotonicity of $ the second term is nonnegative. Hence we should choose C in such a way that the coefficient 
in front of if in the third term is nonpositive. If C > Um, we have 

(3<^{U -C) + a< -/3z7vn^(oo)e^^(^-^«) - ^ Yl + ^• 
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Thus, if it happens that 

-/3 Zjn°{oo) + a< l3zNn%{oo), (20) 

we indeed take C = Um and the corresponding term is nonpositive. If (pOj) is not true, then our choice is 

C = UM + — \og\- - o\ \ >Um. 

This choice of the constant C imphes that ip{y) = and it yields the lower bound in (jlSp . 

Let us switch to the upper bound. We now take ^p{y) = {^{y) — U{y) — C)'^ , where C is another constant to 
be determined. We note that by Lemma[2l ip G H^{Yp) and <I>(^')(^ is integrable. It yields 



/ |Vv9|2 dy + (3 [ ($(^') - + C)){^-U- C)+ dy 
Jyf Jyf 

+ / {P^{U + C) + a)ip dy = 0. 



By monotonicity the second term is nonnegative. Again we should choose C in such a way that the coefficient 
in front of ip in the third term is nonnegative. If C + Um > 0, we have 

P^{U{y) + C) + a> -/3zin?(oo)e-^i(^+^'") - ^Yl ^j^'ji^) + ^■ 



Thus, if it happens that 



P Y ZjH^jioo) -a< -/3zin?(oo), (21) 



we take C = —Um- If (|2ip is not true, then our choice is 



This choice of the constant C gives the upper bound in (|18p and the Proposition is proved. □ 

By classical regularity theory for elliptic partial differential equations, we easily deduce from Proposition [3] 
that the solution of the Poisson-Boltzmann equation is as smooth as the data are. 

Corollary 4. Suppose 5 E and a e C^{S). Then ^ G C°°(y». 

4 The limit case of small 13 

In this section we investigate the case of small values of /3 which occurs for small pore sizes. Equivalently, in 
the case of very dilute concentrations, we can scale all concentration coefficients n^{oo) by a small parameter 
which is multiplied to /3. We shall prove that the solution ^ = of the Poisson-Boltzmann equation (jlOp (with 
a subscript f3 to indicate that we study the behavior when /3 — )• O"*") is uniformly bounded up to an additive 
constant which may blow up. 

Lemma 5. Let ^'^ be the unique solution of ^0\). There exist a constant C , which does not depend on /3, such 
that 

W^p - M{^p)\\m{YF) < C\\a\\L2(^s)- (22) 
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Proof. We recall that the variational formulation (or the virtual work formulation) corresponding to (jlOp is, for 
any smooth 1-periodic function (p, 

I V(^dy + /3 / <^{^,3)vdy+ [ a^dS = 0, (23) 

JYp JYp Js 

where the nonlinear function $ is defined by (jlip . Let M be the averaging operator defined by M{g) = 
Jy^ g{y) dy. Taking the test function (p = — 7W(^'^) in (|23l) we get 

/ \V^(s\^dy + /3 [ ($(^/3) -^(A^(^'/3))(*/3-A^ (*/?)) 

JYf JYp 

+ [ a{^p-M{^,3)) dS = 0. 
Js 

By monotonicity of <I> the second term is nonnegative and Poincare inequality yields the a priori estimate (j22p . 

□ 

When J^adS 7^ 0, we expect that Aii^^p) blows up as /? tends to zero. This is already indicated by the 
L°°-bounds from Proposition [3l More convincingly, choosing (/? = 1 in (j23p leads to 

/3 / ^^p) dy = - f adS, (24) 

JYp JS 

which shows that 7W($(\I'/3)) blows up, at least. Remark also that, for /? = 0, the corresponding boundary value 
problem ([10]) has no solution. However, in the case J^adS = 0, is bounded as will be proved in Section H?3l 

4.1 Formal asymptotics 

We first obtain by a formal method of asymptotic expansions the boundary value problem, corresponding to the 
limit ;5 — 7- 0. There are 3 possibilities. 

Case 1: J^a < 0. In this case it is the most negative valence zi which matters. Since — A4{^i3) is 
bounded, we look for an asymptotic formula 

^piy) = a^ + My) + oil), (25) 

with a constant a/3 — )■ +00 and ipo a function independent of /3 and of zero mean in Yp. We will see that further 
terms in the asymptotic expansion come in general with fractional powers of /3. 
The equality ([2^ becomes 

AT 



Therefore, we have 



-;S Vzjn°(oo) / e^^^('^''+^«(^)+°(i» dy = - [ a dS. 
~i JYp Js 



zin?(oo)/3e-^i'''3 / e^i(^o(^')+°(^)) dy 
Jyp 



N 

\jn°(oo)e-"^'^^ 
and, since a/3 — +00, at the leading order (j26|) reads 



+/3 Vzjn°(oo)e-"^'^^ / e-^^(^°(^)+°(i)) = / a dS, (26) 
,_o Jyp Js 



From (1271) we deduce 



zin?(oo)/3e-^i'^^* / g-^i(V'o(s/)+o(i)) dy = [ a dS. (27) 
Jyp Js 



a/3 = i^ + Co, (28) 

Zl 
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where Co is a constant which may depend on /3 but is bounded as a function of /3 — )• 0. This decomposition 
allows us to eliminate the singular part in the expansion. We thus get the following nonlinear equation for 

V'o(y) 



a dS 



e-^iMy) dy 



in Yp, 



Vipo ■ n = —a on S, 



(29) 



■ipo is 1 — periodic and / ipo dy = 0. 

Iyp 



In our approximate formula for we have neglected terms of order 0(/3^~^2/^i)_ jj^ the canonical case of 2 
opposite valencies {N = 2, zi = —Z2), the neglected term is of order 0{fi'^). 

Equation (p9|) does not contain /3. Rather than using Voi it is more practical to use ^^{y) = ipo{y) + Cq. 
Then we have 

^0{y) = ^-^+My) + O{(3'/''), (30) 

^1 



and (fQ is the solution to the boundary value problem 



-A(^o(y) - zin?(oo)e-^i^«(J^) = in Yp, 
Vifo ■ n = —a on S, 
ifQ is 1 — periodic. 



(31) 



We note that by testing ([3T]) by a constant and integrating, we get 



Zl 71-5 (00 



e dy 



Yf 



a dS. 



Consequently, (pQ solves (j29p except that it is not of mean zero. We have the following simple result. 

Lemma 6. Assume that a is a smooth hounded function such that J^a dS < 0. Then problem i31\) has a unique 
solution ipQ G H^{Yp) such that 

are absolutely integrable. 

Proof. The corresponding functional, to be minimized, is 

M^) = \ I iVv^P dy + n?(oo) / e-^i'^ dy + f dS. 



It is strictly convex and the condition a dS < insures the coercivity. The rest of the proof follows that of 
Lemma [2j □ 



Case 2: Jgcr > 0. In this case it is the largest positive valence zn which matters. At the leading order, ([24l 
reads 



/3z7vn^(oo)e-^^'^'^ / e-^^(V'o(y)+o(i))) 



Yf 



a dS. 



(32) 



For the same asymptotic expansion (j25p . equation ()32p allows us to compute the singular behavior 0^3 and we 

get the following equation for the zero-mean perturbation Tpo, i/jq dy = 

Jyf 



a dS- 



-ZNlpO 



Iyf fi-^Ni>o{x) dy 

Vlpo • n = —a on S, 

ipo is 1 — periodic and jy^ ipo dy = 0. 



in Yf, 



(33) 
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By the same reasoning as in the first case, we deduce 

^/3(2/) = ^ + eo(y)+0(/3V--), (34) 

where is the solution of 

-ACo(y) - Z7vn?^(oo)e-^^««(2') = in Yp, 

V^o • n = -a on 5, (35) 
^0 is 1 — periodic. 

By testing (|35]) with a constant and integrating, we get 

ZNn%{oo) [ e-^^«o(^) dy = [ a dS. 
Jyf Js 

Consequently, solves (j33p except that it is not of zero average. We have the following simple result. 

Lemma 7. Assume that a is a smooth bounded function such that J^a dS > 0. Then problem i35\) has a unique 
solution ipo G H^{Yp) such that 

are absolutely integrable. 

Case 3: jgcr = 0. In this case the problem corresponding to /? = has a solution and the analysis is much 
simpler. 

The following limit problem 



f -A^o(y) = in Yf, 

(36) 



V^'o • n = —a on S, 



^0 is 1 — periodic and / <^(^o) dy = 0, 

Jyf 



has a unique solution since the function <I> is monotone. Then we have 

^^{y) = ^o{y) + 0{(3). (37) 

4.2 Rigorous perturbation results when J^a dS ^ 

We focus on the case J^a dS < 0: the opposite one, fgcr dS > 0, is completely analogous. Motivated by the 
discussion leading to ([5U|) . we look for in the form 

*/?(?/) = ^ + ¥'/3(y), (38) 

Zl 



where ipp is the solution of 



-Aifpiy) - zin?(cx))e-^i^^(^) + ^(v?/?) = in Yp, 

Vipjs ■ n = —a on S, (39) 
ifis is 1 — periodic, 



with 



N 

$(g) = - ^ Zjn°(oo)/?i-^^/"ie-"^^. (40) 
i=2 

We start with a uniform if -^-estimate for 99^. 

Lemma 8. Let a be a smooth bounded function such that J^a dS < 0. Then, for small enough /3, the solution 
ipp of (E^j satisfies the estimate 

\Wp\\m{YF)<C^ (41) 

where C is independent of p. 
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Proof. The variational formulation of problem (j39p reads, for any smooth 1-periodic function if, 

[ V^^-V^dy+ [ (l>(v3/3) - zin?(oo)e^^i'^'9) V3dy + / aip dS = 0. (42) 

JYp JYp ^ ^ J S 

In ([32]) we take ip = ip/s = ip^ + 99^ and we get 

\V^P^dy + iz^f4{oo) [ \^+\^ dy + (3^^-^)Mi^p 

Yp JYp ^1 



+ / aip+dS+ [ \V^~\^dy + m^^)M{^^) 

JS JYp Zi 

+ j a{^-^-M{^-^))dS+ (^j adS^ M{^]^) < 0. 



(43) 



Indeed, 



and 



/ \Vip^\^dy= [ |V^+|2dy+ / \V^Z\^dy, (44) 

JYp JYp JYp 

(45) 



j aippdS = j a(f^dS + j aiip'^ - M{ip^)) dS + (^j adS^M{ip 
Furthermore, both functions ^{g) and g — ?• — zie~^^^ are monotone and 



$(0) - zin?(oo) = /3$(i^). 
Thus, we deduce 

(l>(^/3) - zin?(oo)e-^i'^'' - $(0) + zin?(oo)) > 0. (46) 
However, we use a further argument of strict monotonicity for — zie~^^^, namely 

(-zin°(oo)e-^i'^'3 + ^in°(oo)) (^+ = (^-^in?(oo)e"^^^^ + zin?(oo)) 



> 



(zi)2n?(oo)(/.+ , (47) 



because (—2:16 ^1^)' = (21)^6 > (zi)^ for 5^ > 0. Equalities (f44]) . (fiSl) . together with the lower bounds (f46]) . 
(jl7|) . applied to the variational formulation (fl2|) . yield the desired inequality (fl3l) . We recall that 

lim /3$(i^) = -zinO(oo) > 0, 

so that, for sufficiently small /3 > 0, /3$(i^^) is a positive bounded constant. Further, the product (/^ fi dS)M.{ip'^] 
is nonnegative. Therefore it suffices to apply Poincare inequality and ()4ip follows. □ 

Next we need a uniform L°°-bound for (^^, as /3 goes to 0. (Recall that the L°°-bounds of Proposition [3] are 
not uniform with respect to /3.) 

Proposition 9. For sufficiently small /3 > 0, we have the hounds 



U{y) -Um- — logmax{l, ^- — -} > ipp{y) > 



U{y) -Um-- logmin{l, (48) 

zi zin{{oo) 



where U is the solution of the Neumann problem |j7[j. 
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Proof. We start with the variational formulation for Lpp — U which reads, for any smooth 1-periodic function 99, 



V{ipp -U) -VLpdy - zi'di{oo) I e~''^'^^(pdy + <^{(fp)(fdy 

Yp JYp JYp 



+W / (pdy = 0. (49) 
JYp 

We take ip{y) = ((/9^(y) — ?7(y)+C)~, where C is a constant to be determined. By virtue of LemmaO G H^{Yf) 
and are integrable. Since the function g — )• —Zje~^^^ and g — )• $((7) are monotone, we deduce from (j49|) 

j \V^\'^dy + j (^-zin?(oo)e"^i(^^^) + $(C/-C) + a)v?dy <0. 

Hence we want to choose C such that the expression in front of ip in the second integral is nonpositive. We have 

-zin?(oo)e-"i(^(^)-^) + ^(U{y) -C) + a< 

Now if 

< 2171,5(00) < 0, 

we take C = Um and the left hand side of (j50p is nonpositive for /3 sufficiently small because the sum in (jSOp is 
small. If not, then our choice is 

C >Um + — log ( ^— ) > Um- 

This choice of the constant C implies that 92(2/) = 0, for small enough /3, and yields the lower bound in (j48p . 
For the upper bound we take (/?(y) = {^piy) — U{y) — C)~^, where C is a constant to be determined. It yields 

j \V^\'^dy + j (^-zin?(oo)e-^i(^+^) + $([/ + C) + a) v?dy < 0. 

Hence we should choose C such that the expression in front of 99 in the second integral is nonnegative. We have 

-2in?(oo)e-"^(^(j^^+^) +^{U + C) + a> (51) 
-zinO(oo)e-^i(^'"+^) - ^ /3i-^^/^izjnO(oo)e-^^(^'"+^) + a. 

Now if 

zin5(oo) < a < 0, 

we choose C = —Um and, for sufficiently small /3, the right hand side of ([5T]l is positive because the sum over 
j G is small and the expression in front of if in the second integral is nonnegative. Otherwise, we choose 

C > -Um - — log ( > -Um 



z\ \z\rv{\po) 

and, again, for sufficiently small /?, the right hand side of ([5T]) is positive, which implies the upper bound in 
is]). □ 

As an immediate consequence of Proposition [9l taking the limit as /3 goes to 0, we obtain the following 
corollary. 

Corollary 10. Let tpo be the solution of i31\}. It satisfies the -estimate 

U{y) - Um - — logmax{l, ^ — -} > ipo{y) > 

z\ z\rv{\po) 

Uiy) - Um - -logmin{l, ^}. (52) 

zi zin\[oo) 
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Theorem 11. We have 

\\<f0-Vo\\cHY^))^C(3'-^'/'\ (53) 
for every positive integer k. Furthermore, let ipi be the solution for 

-A(^i + (zi)2n?(oo)e-^i'^Vi = Z2nl{oo)e~'^'f''' in Yp, 

Vipi • n = on S, (54) 
(fi is 1 — periodic. 

Then, for every positive integer k, we have 

\Wp -^0- P'-^'^^'McHY^)) < (55) 

where < q = min (l — z^/zi , 2(1 — Z2/Z1)). 

Proof. First we observe that ipi3 — (po satisfies the variational equation 



- ^0) -V^dy- zmlioo) / {e-^^^^ - e-''^")^dy = 

Yp JYf 

— / ^{(pp){pdy for all smooth 1-periodic ip. (56) 

JYp 

Now we take (p = (pj^ — ipQ as test function, use the strict monotonicity of the function g — )• —Zje~^^^, the 
L°°-bounds (gSl) and 1^ to conclude that 

y^-MmiYp)<C(3'-''/'\ (57) 

Next we write the equation for (pj^ — ipQ as 

-A{(pi3 - ipo) + {(pfj - (po) = {<Pf3 - <po) 

+zin?(oo)(e--i^f -e-^i'^«) -l>((/?^) in Yp, 
V{ipi3 — cpo) • n = on 5", 
(99^ — ipo) is 1 — periodic. 

Using the estimate (j57p . we get the error estimate of the same order. After bootstrapping we obtain the 
required error estimate (|53p. 

Eventually, we write the equation for ^p — ^PQ — ii^~^'^l'^^'p>\ and repeating the above procedure yields (j55p . □ 

Remark 12. In the frequently considered case of two ions of opposite unit charge (N = 2, —zi = Z2 = 1), 
normalizing the coefficients n^{oo) = "^2(00) = 1, we have 

$(5) = 2sinhg , = + + + • • • 

and, in the case f^a dS < 0, the equations for the functions ipj read 

-A(^o + e'^" =0, 

and we have 

\Wp -^0- /3Vi - /3V2|lc^(y^)) < (59) 

for every positive integer k. 

The case J^o" > is analogous and it is enough to repeat the above strategy with zi replaced by zn- 
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4.3 Rigorous perturbation results in the case a dS = 

Here the proofs are much simpler than in the previous subsection. We just state the results. Again, the starting 
point are the uniform //^-estimate for 

Lemma 13. Let a be a smooth function such that J^a dS = 0. Then the solution of problem ilO\} satisfies 
the uniform estimates 

1 I z -n^ioo) \ 

U(y) - f/™ - - logma. [l, - j > > 

U{y) -Um log max 1, - V ' i\ \ , (60) 

ZN y Z7vn'^(oo) j 

and 

\\^p\\m(Yp) < C, (61) 

where C is independent of p. 

Proof. The L°°-bound (j60p is a direct consequence of Proposition [3l Note that (j60p is uniform with respect to 
/?. To obtain (|6ip we take the test function ip = in the variational formulation (j23p 



/ iV'^fsfdy + p / ^{-^ 13)"^ 13 dy + / a'^fsdS = 0. 

Jyp Jyp Js 



'Yf JYp 

Since $ is monotone and satisfies $(0) = 0, we have ^{'^ p)'^ p > 0, while the assumption J^cr dS = implies 
that 

/ a-<l>pdS= / a{^p- M{^p))dS <C\\a\\L2^s)\\^^p\\L^{Yp) 
Js Js 

by virtue of Poincare-Wirtinger inequality. We thus deduce 

||V^'/3||L2(yp,) < C||fT||i2(5), 

which, together with (j60p . implies (j6ip . □ 



In subsection 14.11 we already introduced the limit problem ()36p for ^'o = Ivca.p^Q'^ p. We can also define a 
corrector as the unique solution of 

' -A^'i(y) = -^^o) m Yf, 

^ V^i • n = on 5, ^g2) 

^'i is 1 - periodic and / $'(^o)*i dy = 0. 

JYf 

There exists a solution of ([62]) because /y^ *I*(^o) dy = as required by the definition of (j36l) . As an obvious 
consequence of Lemma [13] we get the following error estimate. 

Theorem 14. Let be the solution of [3^j and ^'i that of (6^1 . Then we have 

11^/3 - ^o|lc^(y») <C/3, ||^';3-^o-/3^'i||cfc(y^)) <C/32, (63) 

for every positive integer k. 

The proof of Theorem [T3] follows the lines of the proof of Theorem [TTl 
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5 Large (5 limit 



We now investigate the asymptotic behavior oi p when the /3 parameter goes to +00. In view of its definition 
dni), a large value of /3 corresponds either to a large pore size L or to a small Debye length Ad, but also to a large 
common value of the concentrations n^(oo). A similar asymptotic analysis has been performed in [2], [16] in one 
space dimension. In higher space dimension our main tool to obtain the behavior near the solid boundaries is 
the multidimensional boundary layer technique introduced by Vishik and Lyusternik |22j . 



5.1 Formal asymptotics 

In the Poisson-Boltzmann system (llOp the parameter (3 appears in the partial differential equation but not in 
the Neumann boundary condition. This indicates the presence of boundary layers in the asymptotic analysis, 
the thickness of which shall be of the order of 0(l//3). The usual technique to handle this situation is that 
of matched asymptotic expansion. We first consider an outer expansion of the solution ^'^ in Yp, away from 
the boundary S. In a second step we shall construct an inner expansion of ^ p in the vicinity of S, which is 
equivalently a boundary layer. 

We begin with the outer expansion for ^ p which reads 



= ^oo + fl*l,oo + ^2,oo + . . . . 



1 T 1 

After plugging this ansatz in the Poisson-Boltzmann equation (llOp we get 



--^A^^ + $(^^) = $(^oo) + ^(^'(^00)^1,00 -A1'oo)+---=0 in Yf 

which implies, at the zero order, that <I>(^'oo) = 0. The electroneutrality condition (jlSp tells us that is the 
unique root of the monotone function Therefore we deduce = in Yp- In other words we have 

^,3(2/) = ^(^) ™ ^way from the boundary S. (64) 

In fact we will check rigorously in the next Subsection that this order of magnitude holds for the L^— norm of 

We now turn to the inner expansion of ^ i.e., its behavior close to S", which is more complicated. We study 
it locally near a point G S, using the geometrical setting introduced in Subsection 13. 11 We consider a tubular 
neighborhood Yp of S with /_f small but much bigger than /3~^/^. Locally, in a neighborhood AA(yo) of we 
make the change of variables y ^ q = {q' ,qd)-, as defined in Subsection 13.11 which satisfies |Vygrf| = 1 in M{yQ) 
and n • Vyqd = 1 on M{yo) n S. The Jacobian J (corresponding to the volume differential change dy = Jdq) is 
defined by 

J = detf^) , (65) 

and the metric matrix (corresponding to the transformation y ^ q) 



J dyj 



l<k,r<d 



which satisfies 

Kd,d = 1, Kk,d = for 1 < A; < d - 1. 

Notice that the coordinates q = q(y) are introduced in such a way that the level sets {qd = const} and the 
normal lines {q' = (Ci, . . . , Cd-i} are orthogonal (that is the corresponding tangential hyperplanes and lines are 
orthogonal). Since Vyqdiy) gives the direction of the normal line and Vyq^, k = 1, 2, . . . , d — 1, form a basis in 
the tangential hyperplane, we have K^^d = ^yqk • Vydist(y, S) = for k ^ d. 
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Differential operators in new coordinates transform as follows: 

_5_ ^ dqkd_ . ^ ^ ^ 



(9^ 1 

V-^ = -divg(JKV,) mY^nJ\f{yo); (67) 



n • V„ = — — — on b. 

dqd 



Applying this change of variables to the variational formulation (|23|) of the Poisson-Boltzmann system (|10p yields 
the following equation in the new coordinates 

-diYg{JKVq^p)+PJ<^{^l3) = 0. (69) 

Dividing ([69|) by J yields that the partial differential equation in Yp n A/'(yo) and the boundary condition on 
S n M{yo) of (fTOl) transform into 



/3 



„ + /3<l>(^«) + lower order derivatives in 

+ second order differential operator in q' = in Yp nM{yo), (70) 

^ = a on 5nAA(2/o). (71) 
dqd 

As usual in the method of matched asymptotic expansions, problem (I70p - ()7ip serves to construct the inner 
expansion. Since we expect the thickness of the boundary layer to be of order 0(1/13), we search for the inner 
expansion of the form 

00 

^p{q',qd) = r^'^Y.P'^'^'^M'^P^'^'id). (72) 

j=0 

Expanding the zero order term <I>(^'^) in Taylor series and taking into account the bulk electroneutrality condition 
(fT5]l . $(0) = 0, we obtain 

N 



k=l 

N N 

= r'l^ ^i"°(°o)^o + ^'^°(°«) f^i - o ^^(^o)') + . . . (73) 



k=l 

N N 
k=l k=l 

Introducing ^ = fi^/'^qd, substituting ([72|) and ([73|) in ([70]) - ([7T]) and collecting power-like terms in the resulting 
equations, after straightforward rearrangements we arrive at the following problem for the main term of the 
expansion: 

(fi f ^ \ 

;^^o(9',0 - iTzlnlioo)] ^o{q',0 = for ^ > 0; (74) 
^o{q',0)=<r{q') for ^ = 0. (75) 



d_ 



Problem (j74p - (j75p is a second-order ordinary differential equation on the positive half-line. After matching with 
the outer solution ^ p = 0(— ), we impose additionally that ^0(9') +00) = 0. The exact solution of ([7i|) -([75 |) is 
thus 

^o{q\i) = ^^exp{-ev^}, (76) 
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with $'(0) = zf.rrl{oo). Back to the original variables, we get 



^/3||Li(y^) < (80) 



^^iy) = -^0^eM-d{y)vWW)} + Oi^), (77) 

where d{y) is the distance between y and S. This asymptotic expansion (|77p will be justified rigorously in the 
next subsection. 

5.2 Rigorous error estimate 

We start with two useful simple inequalities for the nonlinearity <I>. 

Lemma 15. Let <I> be given by ill]) , i.e., <l>(x) = — Zkn^{oo)e~^''^ . There exist positive constants H,Co,Ck 

such that, Vx G R, 

<l>{x)sign{x) > H^\x\, (78) 
$'(x) > Co + Ck\x\''~^, k>2. (79) 

Proof. To prove (f79|) we note that 

«''(x) > (zi)2n?(oo)e-"i" + (z7v)'n^(oo)e-^^" 

with zi < < ^AT, which implies the desired result. Then ()78p follows from a Taylor expansion of ^{x) at and 
the bulk electroneutrality condition (llSp . $(0) = 0. □ 

We now prove a priori estimates which improve that of Lemma [5l 
Lemma 16. Let ^'/j 6e i/ie unique solution of U0\) . There exists a positive constant C such that, V/3 > 1, 

C 

||*/3||L^(y,) < Cr'/^2fc)^ ^>2, (81) 
||^/3||Hi(y,)<C/3-i/^ (82) 
Proof. First, in the variational formulation (j23p we use the test function ip = p. It yields 

/ \V^p\^dy + /3 [ ^^p)^pdy = - I a^pdS. (83) 
JYf Jyf Js 

Since ;5 > 1 and $(a;)x > ff^ |xp because of ([78]) . we deduce from ([83]) that ||^/3||j|/i(y^) < C (this estimate will 
be improved later). Second, in the variational formulation (j23p we use a test function which is a (monotone) 
regularization of the sign of ^'^3. The first term in (j23p is non- negative and the right hand side is bounded thanks 
to the previous estimate. Thus, after passing to the regularization parameter limit, we get 

/3 / ^{^p)sign{^p)dy<C, 

JYp 

and after applying inequality ()78p we get (j80p . Next, we consider again (j83p where the nonlinear term is bounded 
from below using inequality (j79p 

/3 / ^{^>p)^pdy>p(co [ \^p\''dy + Ckl dy\ . (84) 

Furthermore, using a trace inequality [TI], we get 

< C {m^^pWh^y^^ + (/35)^'/'ll^/3|li/iV,)) , (85) 

where we used Young's inequality ah < a^/A + 3b^/^/4 for a = (/35)i/^||^'/3||J^'^2^(y^) and b = {psy^/^W'i^ (sWh^Yp)- 
For /3 > 1 and 5 > small enough, (j8ip - (j82p is a direct consequence of (|85p. □ 
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Since S is compact there exist finitely many points G S and neighborhoods M{'ifl), 1 < i < M, such that 
the open sets Wi = YFr\M{y^i) cover 5, i.e., S C ufi^TFi. Take Wq CC Yf so that Yp C ufio^^i and let {Ci}»=o 
be an associated partition of unity. Here Yp and S are considered as subsets of the unit torus T*^, so the functions 
Ciiy) are 1-periodic. 

Proposition 17. In each set q{Wi) (the image ofWi by the map y ^ q) define a boundary layer function 

^f{q) = r'^Hoiq',P'^'qd) = exp{-grfvW(0)}, (86) 

where is defined by (TS^ . For any smooth test function ip, such that tp = on dq{Wi) \ q{S Ci Wi), it satisfies 

[ K'Vq^f ■ Vgp fdq + p[ ^{^f)p f dq 

Jq{W,) Jq{W,) 



aipJl + \Vgn^\^dq' = / Rlpdq, (87) 

where ,P is the Jacobian of the map y ^ q, defined by i65\) . is the metric matrix defined by / 1 6' 6]) and 

C 

Proof. By direct calculations, using the explicit formula (j86p and taking into account that i^^ d = ^ ^ k < 
d-1. ' □ 

Considering the boundary layers tpf as functions of y now, we immediately obtain the following corollary. 
Corollary 18. The boundary layers i^Wqiy)) satisfy, for any p G H^{Wi), 

I yy{Cii^i')-yy^dy + P [ ^{Qiijf)^dy+ f da^pdS = 
JWi JWi JsnWi 

BJ-pif dy, 



Wi 



where the redefined reminders i2g satisfy 



C 

l|-R^llL°°(Wi) < C, 11^^11x1(14^^) < (90) 

Due to the geometric assumptions from Section 13. H for sufficiently small, the tubular neighborhood Yp = 
{y^Yp : dist(y, 5) < ^} satisfies C uf£iWi. 

An arbitrary smooth function /, defined on Yp is then written as f = Ylf^i Ci/ + Co/- The boundary between 
Yp and Yp \ Yp is and we extend smoothly / from Yp into Yp \ Yp. In Yp \ Yp, /, together with its 
derivatives, is exponentially small with respect to 1//3. 

In the calculations which follow we replace / by the above extension of Cif from Yp to Yp. The error 
is exponentially small in and we ignore it. 

Collecting together the boundary layers with the associated partition of unity, we define 

M 

V''' = EC^V'f (91) 

i=l 

and deduce the following result. 

Proposition 19. For any (p G H^iYp), we have 

I - V'') • V99 + /3 / U{^p) - ^>(^^')) pdy= I Rpif dy, (92) 

JYf JYp ^ ' JYp 

where the global reminder Rp satisfies 

Q 

\\Rp\\L^{yF) < C, ||i?/3||j^i(y^) < (93) 
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Proof. We obtain (j92p by subtracting the variational formulations (j89p of the boundary layers ipf from the 
variational formulation (I23p of We use the fact that 



i=l 1=1 

since = 0{/3~^/^) and $(0) = 0. □ 
Finally we obtain the main result of this section which is a rigorous justification of ()77p (recall that ipi\q) = 

Theorem 20. Let 6e the unique solution of I110\) and ip^^ be given by liQl]) . There exists a positive constant 
C such that, V/3 > 1, 

11^/3 -V'''llLi(y,)<^, (94) 



^3/2' 

c 



l^^-V''''llL^(y.)<^, (95) 



11^/3 -V'''llHi(y.)<^. (96) 

Proof. The proof is similar to that of Lemma [TBI First, we test (j92p by the regularized sign of - iJj^K After 
passing to the regularization parameter limit and using the second inequality of (j93p . we get 



/3 



^ (cI.(M/^) - ^^p'')^ sign (m/^ - V/^) dy < 



Since dZH) implies that $'(x) > C > 0, we deduce (pil) . 
Next we test ([92]) by ^-^^ - V*'- It yields 



< ||i?/3||L-(y,)||^/3-V''''||Li{y,) < (97) 

For /3 > 1, (j95|) - (j96|) is a direct consequence of (j97|) . □ 

Theorem [20] justifies the approximation (I77p by providing error estimates in integral norms. The next result 
gives pointwise estimates for the same asymptotic approximation. 

Lemma 21. There exist positive constants /3o, Ci and C2 such that, for all /3 > /3o and for all y £ Yp, the 
following estimates hold: 

|vl'^(y)| < -^exp{ -C2v^dist(y,5)}, (98) 
|V>I/^(y)| <Ciexp{-C2V^dist(2/,5)}, (99) 
\^(s{y)-i^'"{y)\ < y exp{ -C2v^dist(y,5)}. (100) 
Proof. Introduce a function of s G M 

' $(s)/s for s / 0, 

Pi^) = { ^ 2 0, ^ 

yZjUjlooj for s = 0. 
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From ()78p in Lemma \T5\ we deduce p{s) > > for all s G M. It also readily follows from the definition of p 
that p is a continuous function of s. 

Recall that Yp and S are subsets of the unit torus T'^. Therefore, 1-periodic boundary conditions are implicit 
for all boundary value problems below. For the sake of brevity we do not indicate them. Introducing a function 
Bp{y) = p{y)) (which is continuous and satisfies Bp[y) > H"^), the Poisson-Boltzmann equation (|1U|) can be 
rewritten as 

r -A^';3 + /35/3(y)^/3 = 0, inyj., 
[ • n = —a, on S. 

Denote S = ||c||l°o(5). Then, by the maximum principle, 

< "^fi in Yf, 

where \I'^ is the unique solution of 

r -A*+ + /3i/2^+ = 0, in Yf, 

\ V^'/3 • n = S, on 5. 



Thus, it suffices to show that 

<;^exp{-C2V^dist(y,5)}, for all y G y^. (101) 

In order to prove (jlOip . we are going to construct a so-called barrier function. For any y £ S denote by R{y) 
the radius of curvature of S at y. Under our standing assumption on the smoothness of S we have R{y) > Rq > 0, 
y y € S. For Rq small enough, each of the neighborhoods M{yf), 1 < i < M, covering S contains a ball of center 
y^ and radius {Rq/2). In this ball, we rewrite the Laplace operator in terms of the new coordinates q = q{y), 
introduced in Subsection 13.1^ 

with regular bounded functions Qij{q) and Q^j{q) defined in terms of Kij and J. Setting 

G{s) 

and 



«2 



4^, otherwise. 



after straightforward computations we obtain 



VC/- n 



dU 
dqd 



= S. (102) 

gd=o 



Notice also that U G H^{Yf) because VU{y) = if dist(y, S) = Rq/2. 

Substituting U in the equation yields, for all y G Yp such that q^ < Ro/2, 

Clearly, for all sufficiently large /3 the above right hand side is positive. Therefore, for (3 > (3q, 

-AU + /SH'^U > in Yf. 
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Combining this relation with (jl02p . and using the maximum principle, we conclude that ^'^ < U in Yp. Since 
dist(?/, S) in Yp is bounded by some constant C3 and Q{s) > C4S for any s G (0, C3), this implies the first desired 
estimate ([98]) . 

Estimate (fMj) follows from (f98]l thanks to the standard elliptic estimates. Indeed, in the rescaled coordinates 
z = y/Py equation PU]) reads 

-A,^p = -H^p) in ./pYp 
^z^p ■ riz = \/liS. 
By ([98]) for any ball Bz^^i = {z : \z — zq\ < 1} with yo £ VT^Yp we have 

VP vP 

z G -6^0,1 1^ V^Yp. Considering our regularity assumptions on S and a, by the local elliptic estimates for Poisson 
equation including those near the boundary, we obtain 

|V.*^(z)| < ^e-^^^ii^t^^-^'v^^), z G Bz,,i/2 n y^Yp. 



In the coordinates y this yields the desired estimate ([99p . 

Estimate (llOOp can be obtained by means of similar arguments as in the proof of ([98p . Here we just outline 
the proof and leave the details to the reader. Prom the definition of and tjj^'' and estimate (|98p it readily 
follows that the difference Va = ^'/j — satisfies in Yp, jj, = the following problem: 

-AVp{y) + P^\G)Vp{y) = 9i{y) in Y^, 

with 

\gi{y)\ < c^e--'^^^^^\ \92{V)\ < -^cie-^^v^^«/2 < Icie-'^^v^^"/^. 

here the constants C2 > and C3 > do not depend on f3. Setting Vj3 = ^e"^''^^'''*^^^ and choosing large enough 
C > and small enough Hi > 0, we obtain that for all sufficiently large /3 it holds 



2 



Therefore, V/^ < Vp in Y^. Similarly, Vp > -V ,3 in Y^. This yields (HUOj) in F^. In \ Y^ follows by the 

maximum principle. □ 



6 Dirichlet boundary condition or C, potential at the boundary 

The previous asymptotic analysis was specific to the Neumann boundary condition (or given charge density a) 
imposed on the pore walls 5. The situations is quite different for Dirichlet boundary condition (or Q potential) 
on S. We briefly investigate this case. We modify the boundary condition in the Poisson-Boltzmann equation 

-A^fi + = Q in Yp, 

= C on S, (103) 
p is 1 — periodic. 

All unknowns and parameters are exactly the same as the ones in Section [3.2[ 
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6.1 The limit case of small /? 

We start by studying the behavior of ^'^ when /3 goes to zero. Performing a formal asymptotic expansion 

^'^ = ^0 + + • • • , 
it is easy to check that the zero-order term is constant 

^o(y) = C, 

while the first-order term is the solution of the linear problem 

-A^i = -^'(C) in Yf, 

^-1 = on S, (104) 
^'i is 1 — periodic. 

It is not difficult to justify this ansatz and to prove the following error estimate. 
Lemma 22. There exists a positive constant C such that 

\\-i'p-C-P^i\\miY^)<cp^. 

6.2 The limit case of large (3: formal asymptotics 

We are now interested in the behavior of for large /?. As in subsection 15.11 we are going to use matching 
asymptotic expansion. Of course, the outer expansion, being independent of the boundary condition, is the same 
and we get again (j64p . namely 

^ pin) = 0(— ) away from the boundary S. (105) 

To the contrary, the behavior close to S differs significantly from Section [5l We study it locally, in the same 
geometrical setting as before. We obtain the same differential operator close to the boundary with a different 
boundary condition 



52^ 



/3 



„ + P^i^ r) + lower order derivatives in 
dqj 

+ second order differential operator in = in Yp riAf{yo), (106) 
^p = Ciqi,...,qd-i) on 5nAA(yo). (107) 

The boundary condition (I107p (given ^ potential) is a much stronger constraint than the previous Neumann 
condition (j7T]) . Consequently, we change the inner asymptotic expansion which, instead of (j72]) . is now of the 
form 

^M^', Qd) = ^odl'^ P^'^ld) + /3"^/'*i,c(g', P^'^qd) + .... (108) 

Introducing ^ = /3^^'^qd, substituting (jlOSp in (I106p - (I107I) and collecting power-like terms in the resulting equa- 
tions, we arrive at the following problem for the leading term of the expansion 

-^^0.C(9'' + H^o,dQ,0) = for C > 0, (109) 

^o,c(9',0) = C(?')- (110) 

After matching with the outer solution "J^^ = 0(/3~^), we impose additionally that ^'o,^((7', +00) = and the 
square integrability of the derivative. 
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N 

Let C{x) = ri^(oo)e~^^^ be the primitive of Then the problem (|109p - (|110p admits the first integral 

- ^(^*o,c)^ + C(^'o,c) = Ci= constant. (Ill) 
As we impose ^'0,(^(9', +00) = and the square integrabihty of the derivative, it fohows that 

N 

Ci = C(0) = ^j(oo) > 0. (112) 
i=i 

Thus, plOD . (fm]l and (fTT2]l give 

' ^o,ck=o = C, 

^^o,c = -2sign(C)^C(^o,c)-C(0). ^^^^^ 

Proposition 23. The Cauchy problem Iill3\) has a unique smooth solution ^'o,( on (0, +00), satisfying U09\) 
and 

l*o,c('?',e)l<lC(9')|e-^^ 

4*oc(?^ol<v^lc(gOl^/^e-v^«/^ ^''^^ 



i(;/iere Cs = min{^^.g^._ zjn^-{oo), X^jgj+ 2|n^(oo)} and Cq = 2 max |<i)(C((7'))l ■ 

Proof. For C = 0, the unique solution is ^o,C = ^^"^ there is nothing to prove. 

For C 7^ 0, the Cauchy problem (jll3p has a unique maximal smooth solution on some interval (0,^). If ^ > 0, 
then the solution is positive, monotone decreasing and it reaches the value at ^ = For < 0, it is the 
opposite situation. But is a critical point of (jll3p and no trajectory can leave or reach that point. So the 
solution cannot be zero for some finite L Therefore, the Cauchy problem ()113p has a unique maximal smooth 
solution on the entire real line (0, +00). 

Next, a simple calculation gives 

4^'o,c(9',Ol'>C.|^'o,c(9',e)l', 

with Cs = z|n^(oo), X]j6j+ z'jn^j{co)}. Consequently, for C > 0, we have 

d 



^*o,c(9',0<C(g')e- 



and we establish the exponential decay of ^'o,( and the first part of (jll4p . For C < everything is analogous. 
The ordinary differential equation (IllSp gives 

|^^o,c(9',e)P<2max|$(C)||^o,c(9',OI 

and we conclude the remaining part of estimate (I114p . □ 

Remark 24. In many situations, for a symmetric electrolyte with ion charges ibQ, the explicit solutions are 
known. A classical reference is the book [21^. For example, in the case —z\ = 1 = Z2 and n5(oo) = 712(00) = 1/2, 
we have the following Gouy- Chapman solution 

T / / 01 l + tanh(C/2)e-« 
1 — tanh[Q / 2)e * 

By direct computation we can check that this solution satisfies the properties established in Proposition \23l In 
the general case ^o,( can be expressed using elliptic functions. Nevertheless, our simple analysis gave us the 
properties of the solution without using its explicit form. 
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6.3 The limit case of large /3: rigorous error estimate 

We start with the study of the boundary layer function ^'o.c;- As in section Owe use the local change of variables 

y ^q- 

Proposition 25. For each open subset Wi = Yp riM{y^), let the boundary layer function be defined by 

i^f{q) = ^o,dQ',QdV^) (115) 

for the given boundary data ( on S D Wi. Then, for any smooth function ip such that ip = Q on d(\{S Pi W.-,), we 
have 

[ • V,<^ J' dq + P f dq= [ R\Vg>ipj' dq 

+ / Ri(pJ'dq= [ Riiprdq+ [ _ a{^dq', (116) 

JciiWi) Jq(W^) Jq{dWi\{SnWi)) 



where J* is the Jacobian of the map q — )• y, defined by I165\) . is the metric matrix defined by ( fggi) . Vg' is 
gradient with respect to q' and 

y^{\\Ri\\L-°(q{Wi)) + h\\\L--{q(dWi\(srm,)))) + ^ C^/p] (117) 

\\R]\\Lr(qiw.)) < C/3(i-i/'^)/2 j = 2,3. (118) 

Proof. By direct calculation, using the estimate (|114p . We note that the higher derivatives of ^'o,^ with respect 
to q' satisfy also the estimate ()114p . □ 

Corollary 26. For all ip G H^{Wi) such that ip = on S H dWi, we have 

[ V,Vf -Vyifdy + P [ a>(^f )v9 dy= [ R{VyV dy 
JWi JWi JWi 



+ / R2(pdy= / R\ipdy+ / _ a\ipdS, (119) 



where redefined reminders a\, Rj, j = 1,2,3, satisfy ^117 )- fTTS\) . 



'J' 

Using the definition of Tp^'' , we see that ipf = ijj^f" on Wi n Wj , when Wi H Wj is nonempty. Now let ip^^ = if^f' 
on Wi. As in Subsection 15.21 we make a smooth extension of ip^^ into Yp \ Yp. In Yp \ Yp, together with its 
derivatives, is exponentially small with respect to 1//3. We simply ignore the exponentially small terms in our 
estimates. We have 

Proposition 27. For any (p> € H^{Yp) such that (p = on S , we have 

[ V(^';3 - V''') • + /3 / m^^)-^^p''))ipdy= [ R^V^ dy 

JYp JYp JYp 

M 

+ [ R}ipdy= [ R^ipdy + y^f _<T),^dS, (120) 
JYp JYp ~l JdWAisnw,) 

where the global reminders i?^, j = 1, 2, 3 and o"^, z = 1, . . . , M satisfy 

V^i\\R}\\L-{Yp) + \Wh\\L-{dWi\{SnW,))) + 

\\R^/,\\LriYp)<Cf3^^-^/'-^/\ VrG [l,+oo], j = 2,3, i = l,...,M. (121) 
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Theorem 28. Let he given by [im) and by HlSl) . Th en we have the following behavior for large f3: 

11^/3 -V'''llLi(y^)<^, (122) 

\\^p-i^'%L^^YF)<^, (123) 

11^/3 -V'''llHi{y^)<C. (124) 

Proof. First we test the variant of (jl20p involving boundary terms by the regularized sign of {^p - V^')- After 
passing to the regularization parameter limit, we get 



/3 / {^^p) - ^{ij'')) sign {^p - ^'') dy < C, 
JYf 

and after applying a slight generalization of the inequality (j78p we get (|122p . 

Next we use the variant of (|120p involving only volume terms and test (|12Up hy gp = p — ip^K It yields 



- ^"Or dy + l3 m'^p) - HriM'^p - ri dy < 

Yf JYf 

LHYf) 

/3min$' , ,,r, 1 ,,o ^ , 

< ^^^MhiV,) + 2 + C. (125) 

For /3> f3o, (fT23]l - (fT2i]l is a direct consequence of (fT25D . □ 

Remark 29. By the same arguments as in Section [5l one can obtain, in addition to the inequalities of Theo- 
rem\2Ei pointwise estimates for p and for the discrepancy ^p — 'il)^K In the case of Dirichlet boundary condition 
these estimates read, for /3 > 1, 

\^p{y)\<Clew{-C2^/pdist{y,S)}, 
|V*;3(y)| < Civ^expj -C2V^dist(y,S)}, 

\^p{y)-i^''iy)\ < ^exp{ -C2V^dist(y,5)}, 



where Ci and C2 are positive constants, independent of /3. 
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